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We show that the quantum phase transition of the Tavis-Cummings model can be realised in 
a linear ion trap of the kind proposed for quantum computation. The Tavis-Cummings model 
describes the interaction between a bosonic degree of freedom and a collective spin. In an ion trap, 
the collective spin system is a symmetrised state of the internal electronic states of N ions, while 
the bosonic system is the vibrational degree of freedom of the centre of mass mode for the ions. 
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I. INTRODUCTION 



More than two decades ago, when quantum optics was young, the quantum dynamics of collective spin systems 
interacting with a single bosonic degree of freedom was a major research problem. The model arose as an attempt 
to describe the interaction between a collection of two level atoms and a single mode of the radiation field. Walls 
and co workers [1J] were among the first to realise that such models provided ideal examples of the role of quantum 
fluctuations in the nonlinear interaction between matter and light. Quantum fluctuations were shown to drastically 
change the predictions of scmiclassical theory in such systems. This phenomenon has appeared more recently in the 
discovery of quantum phase transitions in quantum spin glasses J2j and other many body quantum systems. While 
the collective spin models did not directly apply to achievable experiments at the time, they did provide insight that 
subsequently proved important for many other quantum optical experiments including anti-bunching, squeezing ||, 
and cavity QED Q . In this paper we show that the models of a collective spin interacting with one or more bosonic 
modes can now be experimentally realised in modern ion trap systems of the kind proposed for quantum computation 
§,§. An enormous effort has gone into making such systems work at the quantum level, with little interference 
form classical sources of noise, and a number of such experiments exist today. It would thus appear worthwhile to 
reconsider the collective spin models, and the associated quantum many-body effects exhibited by such systems, with 
a view to direct experimental realisation. 

In particular we consider the Tavis-Cummings (TC) model j7j, which can be realised in a linear ion trap of N ions 
with the bosonic degree of freedom appearing as the quantised collective centre-of-mass motion. If each ion is coupled 
£7^" to the- vibrational motion using an identical external (classical) laser detuned to the first red-sideband transition, the 
symmetry is such that the electronic degree of freedom for the ions can be described as a collective spin (N) and 
• i-h , the reversible dynamics is well described by the TC model. The TC model is known to exhibit important nonlinear 
quantum effects including a quantum phase transition || in which the (zero temperature) ground state undergoes a 
morphological change as a parameter is varied and averages of intensive quantities undergo a bifurcation. 
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II. THE TAVIS-CUMMINGS MODEL 



The interaction Hamiltonian for N ions interacting with the centre of mass vibrational mode can be controlled 
by using different kinds of Raman laser pulses. A considerable variety of interactions has already been achieved 
or proposed ||||||. In this paper we consider the first red-sideband transition. The ion is assumed to be in a 
three dimensional anisotropic harmonic potential. Two dimensions are very tightly bound and are neglected. In 
the remaining dimension, an external laser couples the electronic state to the vibrational motion. If the vibrational 
frequency is large enough and the Lamb-Dicke limit || applies the motional sidebands of the absorption of the 
electronic transition can be resolved and a laser detuned below the electronic resonance by one unit of the trap 
frequency can excite the electronic transition by absorbing one vibrational phonon, the additional energy required 
being made up by the laser. We will assume that the laser ( or lasers if a Raman process is used) is sufficeintly strong 
that it can be treated classically. Under these assumptions the Hamiltonian, in the interaction picture, is 

N 

Hi = hnJ2( a ° ( + + « tcrW ) (!) 
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where the coupling constant is £1 = r/flo where rf = E r /(%Muio) is the Lambe-Dicke parameter with E r the recoil 
kinetic energy of the atom, loq is the trap vibrational frequency, and M is the effective mass for the centre-of-mass 
mode. The Lamb-Dicke limit assumes 77 << 1, which is easily achieved in practice. The frequency, f2o is the effective 
Rabi frequency for the electronic transition involved. The raising and lowering operators for each ion are defined 
by a- — \g){e\ and <r+ = |e)(<?|. This sideband transition can be used to efficiently cool the ions to the collective 
centre-of-mass ground state, thus preparing the system in the vibrational ground state 

If the external laser field on each ion is identical (in amplitude and phase) the interaction Hamiltonian is 

= hn(aJ+ + ,L) (2) 

where we have introduced the bosonic annihilation operator a for the centre-of-mass vibrational mode and where we 
have used the definition of the collective spin operators, 

N 
i=l 

where a = X, y, z . Identical laser fields could easily be obtained by splitting a single, stabilised laser into multiple 
beams. The interaction Hamiltonian in Eq (||) specifies the Tavis-Cummings model 0. This model first appeared 
in quantum optics where the bosonic mode is the quantised field in a cavity. However this realisation is difficult 
to achieve experimentally. In contrast the vibrational mode realisation should be readily achieved. The dynamics 
resulting from this Hamiltonian is quite rich. Collective spin models of this kind were considered many decades ago in 
quantum optics PJTcj] . In much of that work however the collective spin underwent an irreversible decay. In the case 
of an ion trap model however we can neglect such decays due to the long lifetimes of the excited states. On the other 
hand heating of the vibrational centre-of-mass mode can induce irreversible dynamics in the system in a manner that 
has not been previously considered, and that is reminiscent of thermal effects in condensed matter physics. 

We are interested in the driven Tavis-Cummings model in which the vibrational mode is subject to a linear forcing 
term which can easily be achieved by a suitable combination of Raman laser pulses, or by appropriate AC voltages 
applied to the trap electrodes H. In this case the Hamiltonian, in the interaction picture, is given by 

H I = hVL{aJ + +(Jj-) + hE(a + a)) (4) 

This may be written in terms of the hermitian canonical oscillator variables X = (a + a^)/v2, Y — —i(a — a^)/v2, 
and the canonical angular momentum variables J x = (J + + J_)/2, J y = —i{J+ — </-)/2, J z = [J+, J_]/2. It takes 
the form 

H = XJ X - YJ y + X X (5) 
with x — E/Q, and we have scaled the Hamiltonian by H — > H/V2FI. This indicates that time is measured in units 

Alsing has shown that the ground state of this system, for weak driving, is a product state in which the bosonic 
mode is squeezed and the electronic states are rotated in the angular momentum space. We provide a direct proof 
of this statement below. However it is first useful to consider the dynamics of the equivalent semiclassical model as 
many of the results in the quantum case can be interpreted in terms of the features of the semiclassical model. 



A. Semiclassical Tavis-Cummings model 

The Tavis-Cummings model represents an interaction between a simple harmonic oscillator and a linear top for 
which there is a classical model which we now define. We choose the classical model so that the equations of motion 
are of the same form as the Heisenberg equations of motion for the quantum model. The classical Hamiltonian is 
defined as 

H = XJ X - YJ y + X EX (6) 

where X, Y are respectively the canonical oscillator position and momentum variables with the canonical Poisson 
bracket {X, Y} = 1, while J7fe are the three components of angular momentum for a classical top with the canonical 
Poisson brackets {J7i,j7fc} = ^Zk^ijkJk- The equations of motion for a canonical coordinate w is given as usual by 
Poisson bracket with the Hamiltonian w — {w, H}. The equations of motion are, 
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X = -J Y (7) 

Y = -Jx - X (8) 
Jx = -YJ Z (9) 

Jy = ~XJ Z (10) 

j z =XJ y +YJ x (11) 

Note that these equations have a conservation law J x + Jy + J 2 = constant. 

We now justify this choice of classical Hamiltonian by noting that the Heisenberg equations of motion for the 
Hamiltonian Eq([|) have the same form as the semiclassical equations of motion with all variables replaced by the 
corresponding operators. We thus see that the semiclassical equations result form taking moments of the Heisenberg 
equations and factorising all product moments. The factorisation assumptions ignores correlations which scale as 1/N 
for the scaled operators J a /N . The conservation law J x + Jy + J z — constant is a reflection of the operator relation 

which in the semiclassical limit indicates that J x + Jy + J z = ^- . 

The classical equations have one nontrivial fixed point at X* = Y* = J* = and J* — — x, J*, — \J N 2 /4 — x 2 - 
However as the conservation law requires that \J X \ < N/2 we see that we must have 

2E 

— < 1 (below threshold) (13) 

which corresponds to an energy of TL = 0. We will refer to this as the below threshold case. As E is increased from zero, 
the fixed point for the angular momentum system rotates about the J v direction eventually reaching the equatorial 
plane at J x = —N/2 at the threshold condition. The oscillator system always has zero amplitude below threshold. If 
we linearise around this fixed point we discover that it is an unstable hyperbolic point with time constant proportional 
to . . Note that this time constant goes to infinity as the fixed point is approached as is typical for a hyperbolic 



fixed point. 

We now consider the above threshold case 

2E 

— — > 1 (above threshold) (14) 

Clearly the value of \J X \ cannot increase above N/2. Indeed there is no fixed point above threshold. However there 
is a special solution curve that continuously joins to the below threshold case for phase curves with H = 0. 

To see this we consider making a canonical transformation by a rotation in both the X — Y plane and in the J x , J y 
plane (see figure [l]). The canonical transformations are 

X = Xcos0 + ysin0 (15) 
Y = F cos 6» -X sin 6» (16) 
Jx = Jx cos 6 -J v sin 9 (17) 
J y = J x cos 9 + Jy sin 6 (18) 

The Hamiltonian then takes the form 

H = X(J x + X cose)-Y(J y - X sme) (19) 
The phase curves with TL = now correspond to cither 

X = ; J y ^ X sin9 (20) 

or 

Y = Q ; J x = - X cos6 (21) 
These phase curves smoothly join the fixed point at threshold if J x = —N/2 which implies 
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(22) 



These solutions are illustrated in figure Note that as E — > oo we have that J x eventually points in the direction 
of — Jy while phase curve in the oscillator phase space points along the Y axis, indicating that for large driving 
the system is essentially a particle in a linear potential which accelerates at constant rate. These results were first 
obtained by Alsing and Carmichael fll^] , 

Y J 




FIG. 1. An illustration of the canonical transformation used in the semiclassical equations above threshold 



III. QUANTUM STATES 

First note that the ground state when there is no driving is \j, — j) z C8> \0) v with a zero eigenvalue. This ground 
state corresponds to the fixed point of the semiclassical model with zero oscillator amplitude and angular momentum 
pointing in the —J z direction. We postulate that as the driving is increased form zero the ground state of the 
Hamiltonian Eq(^) is given by 

\£o)=S(r)R(9)\j,-j) z ®\Q) v (23) 

where \j, — j) z <g> |0)„ corresponds to all ions in the ground state and the vibrational mode in the ground state. The 
operator S(r) is a squeezing operator defined by 



S^(r)aS(r) =iia + va) (24) 



with \i — coshr, v — sinhr. 

The rotation operator R(9) is defined by 



R(0) = e -9{J+-J-) (25) 
and corresponds to a rotation of 29 around the J y axis. Consider now 

Hi\£o) = SR (R^S^HSR) \j, -j) z <8> \Q) V (26) 
If we now transform the Hamiltonian and require that 

RlS^HSR\j,-j) z ®\Q) v = (27) 



we find the following conditions, 



which requires that 



v(l + cos 2(9) = fj,(l - cos 26») (28) 
Qjsin29 = E (29) 



cos26» = e- 2r (30) 



and the ground state energy is taken to be £q — 0. The ground state is thus a product of a squeezed state for the 
vibrational mode and a rotated angular momentum state, rotated about the J y axis. 

The above results are consistent with the semiclassical approximation. The mean amplitude of a squeezed vacuum 
state is zero, corresponding to the semiclassical fixed point at X = Y = while the rotation around the J y axis 
corresponds to the semiclassical fixed point at J x = —\. 

If we continue to increase E above the threshold value the system adiabatically follows a zero energy state, although 
this is no longer a ground state. In fact the canonical transformation used in the semiclassical analysis can be applied 
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to the quantum operator valued Hamiltonian. The result is the same as the semiclassical case, Eq @ with all 
variables replaced with the corresponding operators. The zero energy state then corresponds to the zero energy 
eigenstate of Y cos 6* — X sin.9 with cos6> = NQ./2E. This is of course just a rotated, infinitely squeezed state. The 
electronic state is likewise a angular momentum eigenstate rotated from \ j, —j) in the equatorial plane (orthogonal to 
J z ). Thus above threshold the zero energy eigenstate deforms continuously from the state at threshold. 

Let us summarise these results. For no driving the ground state corresponds to the oscillator in the ground state 
and all ions in the ground state. As the driving is increased, but kept below threshold, this state deforms to a 
squeezed oscillator state while the collective spin system begins to rotate about the J y axis. Note that the mean 
oscillator amplitude (a) remains zero as does the mean of the y-component of the collective spin. As the driving 
increases through the threshold value, this state changes its character so that a non zero value of J y is acquired and 
the oscillator is infinitely squeezed in a direction at an angle cos 9 = NQ/2E to the below threshold squeezing. This 
morphological change of the state as the driving passes the semiclassical critical point is a quantum phase transition. 
The quantum phase transition can be seen in the mean value for J y and J z as shown in figure |^. Below threshold the 
scaled mean values are given by 



and above threshold we have 



where x = 2E/NQ. 
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Vl -x 2 (32) 
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FIG. 2. The scaled monents, (a) 2\(J X )\/N , and (b) 2\(J X )\/N plotted versus the scaled driving strength x = 2E/(NQ). 



What are the experimental manifestations of this transition ? Needless to say no one is ever going to observe 
an infinitely squeezed state in an experiment. So what does happens at 29 = n/2 when the electronic state is 
the J x eigenstate \j — j) x and the vibrational mode appears to be infinitely squeezed ? Is such a state physically 
possible ? Suppose for example we begin in the ground state of the Hamiltonian with no driving (E — 0) which 
is simply \j,—j) x <8 |0)^, and adiabatically increase the driving strength. It would appear that the system would 
then adiabatically evolve into the squeezed vibrational state described above. If we were ever able to reach the case 
29 = 7r/2 we would have reached an infinite energy state for the vibrational mode at a finite driving strength. Clearly 
this is not possible and to understand why it is useful to reconsider the semiclassical dynamics for this model. The 
adiabatic approximation requires that we vary the driving strength on a time scale slower than all other time scales 
in the system. The key time scale for the ground state variation is just the time scale associated with the hyperbolic 
unstable fixed point, (N 2 /A — x 2 ) _1 ^ 2 i which goes to infinity as we approach j& = 1. Thus the adiabatic increase 
of the driving must proceed infinitely slowly, that is it must be switched to the finite value E — in an infinite 
amount of time. This pumps an infinite amount of energy into the system and results in infinite squeezing in the 
centre of mass vibrational mode. Obviously in practice this cannot be achieved so the totally squeezed ground state 
is not possible. However it will still be possible to achieve some squeezing of the vibrational mode at smaller values of 
the driving. This would make an interesting observation for current ion trap experiments even with only a few ions. 
The squeezing of the vibrational mode can be observed using the dynamical method of reference |l3f| 
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In current ion trap experiments, laser cooling techniques allow the centre of mass mode to be prepared in the 
ground state. Unfortunately it does not stay there. Heating due to a variety of sources, including fluctuating linear 
potentials, lead to an irreversible evolution away from the ground state. If such heating is present during the coupling 
of the electronic and vibrational motions, irreversible dynamics will be spread to the collective spin degrees of freedom 
as well. 

As an example we consider what happens if we use the Tavis-Cummings interaction (excitation on first red sideband) 
in the presence of strong heating. Heating of the centre-of-mass mode due to fluctuating liner potentials may be 
described in the interaction picture by the master equation, 

^- = -Xl[aj+ + (fj-,W] + | (V[a] +£>[a t ]) W (35) 

where W is the density operator for the spin and vibrational degrees of freedom and the superoperator T> is defined 
by T>[A]p = 2ApA^ — A^ Ap — pA^ A . The irreversible term corresponds to two point processes in which phonons 
are removed or added from centre of mass mode at the rates j(a'a) and j(aa^} respectively. This does not change 

any first order moments, however it does lead to a diffusion in energy as d ^ t a ^ = 7 . The effect of heating can be 
included in the semiclassical analysis by adding an appropriate stochastic term. In the Ito calculus Hj the effect is 
to add to the equations for X, Y terms of the form 

dX =(...) + VW*(*) (36) 
dY ={...) + ^dW y (t) (37) 

where dWi(t) are independent Wiener processes. If the heating rate is small enough these terms can be neglected. 
However if they are large new steady states can occur in the semiclassical and quantum descriptions which will be 
described in a future publication. 

We would like to thank Howard Carmichael for useful discussions. 
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